AD-A062  261 


UNCLASSIFIED 


STEVENS  INST  OF  TECH  HOPOKEN  N J DAVIDSON  LAB  F/S  8/3 

ON  THE  FEASIBILITY  OF  USING  THERMISTORS  TO  MEASURE  OCEAN  CURREN— ETC (U) 
JAN  62  J W HOLLENBE.R6 


SIT-0L-N-6A2 


5 -78 


, 


II 

II 


>J 

Q_ 

O 

C_5 


< I I I POINT  M \ I ION 
IIOHOKIN.  MW  J i;it>K\ 


Hy drtlffutl  Op+Hr 


4 


6 


DAVIDSON 

LABORATORY 


% ftei/c 


Please- 


i DlSTfcr,,-. 

4Lii‘0l'J  uA 


ON  THE  FEASIBILITY  OF  USING 
THERMISTORS  TO  MEASURE 
OCEAN  CURRENT  VELOCITIES 

by  Joel  W.  Hollenberg 

D D C 

D|aPR  « 1978 

JlkraEtnn 

rJ_^  B 

42  JANUARY  1962 


■ | mm 

I 


Agpwrid  far  public  r*l« 
DfaMbattos  ObIbM 


/ 


\ 


J}N  THE  FEASIBILITY  OF  USING 
THERM I STOBS  TO  JIEASUBE 
(JtEAN  CURRENT  VELOCITIES 

^ s. » _2_ 


by  Joel  * ./Ho  l l 


enberg  j 


Work  Supported  by 
Office  of  Naval  Research 


DAVIDSON  LABORATORY  ^ 
•Stevens  Institute  of  Technology 
Cast le  Point  St  at  ion 
Ilohoken,  New  Jersey 


D D C 

nprarpm  nr? 

UJ  APR  « 1918 


T5E0TTE 


B 


DISTRIBUTION  STATEMENT  A 

Approved  lor  public  release; 
Dtotilbudoe  Unlimited 


TABLE  OF  CONTENTS 


I.  Theoretical  and  Design  Considerations 

II.  Background 

III.  References 


page 


1 

6 

9 


LIST  OF  ILLUSTRATIONS 


Figure 


The  mis  tor'  Probe 


2.  Mechanical  Drawing  of  Thermistor  Velocity-Probe 


3.  Schematic  Diagram  of  Amplifier  Circuit 

4.  Schematic  Diagram  of  Diode  Linearizing  Circuit 

3.  Temperature  Bridge  for  2000  OHM  Thermistor  (Pins  are 
Lettered  for  Sanborn  Input) 


6.  Calibration  for  Oceanographic  Use 

7.  Oceanographic  Data 

8.  Wave  Channel  Data 

9.  Wave  Channel  Data  to  Determine  Sense  Indication 

10.  Wave  Channel  Data  to  Determine  Time  Constant 

11.  Theoretical  Response  of  Fenwall  ffB43  p2  Thermistor 
with  RAT  ~ 106 


1 ACCESSION  for  U 

NTIS 

Whit*  Section 

DOC 

BoN  Section  a 1 

UNANNOUNCED 

□ 

JUSTIFICATION  

BY  .....  — , 

BBIUBUItOK/WMURlUlf  COOES 

1 Dist.  AVAIL  aitf/ot 

SPECIAL 

ftl 

I. 


THEORETICAL  AND  DESIGN  CONSIDERATIONS 
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* The  characteristic  of  thermistors  which  makes  them 
useful  for  measuring  fluid  velocities  is  their  large  negative 
temperature  coefficient  of  resistance,  which  can  be  as  high 
as  -5.8  percent/degree  centigrade  at  room  temperature,  com- 
pared to  0.30  percent/degree  centigrade  for  platinum  (refer- 
ence l) . A thermistor  Immersed  In  a stationary  fluid  can  be 
heated  by  means  of  a current  to  a temperature  higher  than 
its  surroundings.  Thus  It  will  lose  heat  by  steady-state 
convection.  When  a thermistor  Is  exposed  to  a moving  stream 
of  fluid,  it  loses  further  heat  by  forced  convection.  The 
temperature  of  the  thermistor  and  hence  its  resistance  depends 
on  the  fluid  velocity  and  temperature  as  well  as  the  heating 
current.  Assuming  a constant  ambient  fluid  temperature,  a 
relationship  exists  between  the  fluid  velocity  and  the 
voltage  across  the  thermistor.  ^ 

The  velocity  can  be  determined  by  either  maintaining 
the  thermistor  at  constant  temperature  and  adjusting  the 
current,  or  by  heating  the  thermistor  with  a constant  current 
and  measuring  the  voltage  drop  in  the  thermistor  circuit. 

For  fast  response,  constant-temperature  operation  is  superior 
because  it  diminishes  the  effect  of  time  lags  caused  by  the 
thermistor  ( reference  2) . Current  required  for  constant- 
temperature  operation  can  be  adjusted  with  appropriate  feed- 
back amplifiers,  and  the  adjustment  may  be  read  out  as  a 
change  in  bridge  voltage. 

Considering  constant- temperature  operation,  the 
power  supplied  to  the  thermistor  can  be  expressed  as: 


(1) 
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wnere 


V,p  = thermistor'  voltage 
R,p  - thermistor1  resistance 

In  constant- temperature  operation,  R,p  is  constant 
and  Vp,  can  ue  measured.  Tills  power  is  dissipated  in  trie  form 
of  iieat  and  carried  away  by  the  moving  stream.  Therefore, 
an  energy-balance  equation  can  be  written  as 


i(T,r 


(2) 


where 

h = convective  heat  transfer  coefficient 
c 

A = thermistor  surf ace-area 

Tp,  = thermistor  temperature 

Ta  = ambient-fluid  temperature 

If  the  thermistor  temperature  distribution  is 
assumed  to  be  uniform  and  constant,  its  resistance  can  be 
determined  U3ing  data  supplied  by  the  manufacturer.  Because 
the  ambient-fluid  temperature  and  the  surface  area  of  the 
thermistor  are  known,  hf>  only  must  be  determined. 

Considering  the  thermistor  as  a sphere,  an  expression 
involving  hc  for  Reynolds  numbers  less  than  2000  can  be 
written  ( reference  3)  : 

h d 

p^-  Prr“0*3  - 0.97  + 0.b8  Re°  * 5 (3) 
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p - I'luld  density 
V = fluid  velocity 
d = thermistor  diameter 
p = fluid  viscosity 

Equation  5 solved  in  terms  of  nc  yields, 
hQ  = A + B fv  (4) 

and  substituting  equation  4 in  equation  1 yields 

VT  = C(C  + D fv)  (i^AT)]°-5  ( b) 

using  AT  for  the  temperature  difference  between  the  thermistor 
and  the  fluid. 

Pip  and  AT  are  constants  in  a constant  temperature 
mode  of  operation,  and  the  voltage  is  proportional  to  the 
one-fourth  power  of  the  velocity  in  the  Reynolds  number  range 
considered.  Therefore,  Vrp  can  be  increased  at  conatant 
velocity  by  increasing  Ftp  or  AT.  However-,  they  cannot  be 
increased  simultaneously,  because  an  increase  in  Ftp  will  cause 
AT  to  decrease  due  to  the  properties  of  the r-mistors . The 
relationship  between  the  resistance  arid  temperature  of  a 
thermistor  may  be  expressed  as  ( reference  l) : 

-j&p-j-  = exp  /3(ct  - m ) (6) 

' 1 o'  1 o 

where 

R( t)  = resistance  at  temperature  T 
R( T ) = resistance  at  temperature  Tq 
p>  - a constant  of  the  material 
The  best  operating  point  of  a thermistor  can  be 
found-- that  is,  the  largest  voltage  can  be  obtained  for  a 
given  velocity  by  substituting  for  R,.,,  In  equation  b5lts  value 

In  terms  of  Tm  and  T . The  resulting  relationship  between 

T a 

V™  and  T_,  can  be  differentiated  to  find  a maximum  for  VT. 

1 i 

This  sets  the  best  operating  point  of  the  thermistor.  Using 
appropriate  values  for  the  various  quantities  in  equations 
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5 and  6,  a temperature  difference  of  50° F between  the 
thermistor  and  the  water  yields  a maximum  VT. 

For  the  velocity  probe  developed  at  Davidson 
Laboratory,  equation  5 can  be  used  for  velocities  up  to  4.5 
feet  per  second.  Figure  10  shows  the  solution  of  V,p  as  a 
function  of  velocity  obtained  from  these  theoretical  consider- 
ations using  appropriate  values  for  the  various  quantities. 

The  shape  of  the  curve  and  the  order  of  magnitude  of  the 
voltage  are  in  agreement  with  our  previous  experimental 
experience  ( reference  2) . To  extend  measurements  beyond  thl3 
region,  a knowledge  of  heat  transfer  from  spheres  in  the  range 
of  Reynolds  numbers  of  interest  and  an  experimental  program 
of  testing  are  necessary. 

The  use  of  thermistors  as  fluid  velocity  meters  has 
been  tried  by  others  a3  well  as  Davidson  Laboratory  with 
moderate  success  (reference  2).  At  present,  the  main  diffi- 
culties encountered  in  their  use  are: 

(a)  Their  inherent  nonlinearity--the  signal  produced 
by  a thermistor  i3  a function  of  the  velocity 

to  the  one- fourth  power,  depending  on  forced 
convective  heat  transfer  and  flow  regime. 

(b)  Their  sensitivity  to  ambient  temperatures 
(reference  J*)--an  intrinsic  property  of  the 
device. 

(c)  Their  interaction  with  a conducting  1'luid-- 
supplying  electrical  power  to  a thermistor 
in  sea  water  requires  either  insulating  the 
thermistor  to  prevent  electrolysis  or  the  use 
of  alternating  currents. 

( d)  Their  general  insensitivity  to  flow  direction 
(reference  5) --the  cooling  effect  of  the  fluid 
velocity  is  generally  independent  of  the  ther- 
mistor’s orientation. 

(e)  Their  tendency  to  change  In  characteristic 
during  use.  (This  is  a problem  of  hot-wire 
anemometry  also) . 

( f)  Their  physical  fragility. 
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However,  continuing  efforts  to  eliminate  these 
dll  fie ul ties  have  met  with  moderate  success  and  may  be  sum- 
marized as  follows: 

(a)  The  inherent  nonlinearity  in  the  thermistor 
response  can  be  compensated  for  by  appropriate 
electronic  circuitry.  A diode  linearizing 
(reference  5)  circuit  can  be  used  to  rectify 
the  thermistor  response  or  a nonlinear  amplifier 
can  be  used  to  compensate  the  thermistor-  non- 
linearity or  both  may  be  used  together.  In 
addition,  the  thermistor  can  be  operated  over 

a limited  linear-  range. 

(b)  Corrections  for  ambient  temperature  can  be 
achieved  by  using  a second  thermistor  to  measure 
the  fluid  temperature  near  the  velocity  ther- 
mistor, whose  response  at  various  ambient  tem- 
peratures under  a constant  operating  temperature 
can  be  determined  by  calibration. 

(c)  The  thermistor  can  be  shielded  from  a conducting 
fluid  by  coating  the  material  with  a thin  layer 
of  glass  at  a cost  in  response  time  or  sensitiv- 
ity or  both.  Time  constants  less  than  one  second 
and  resolutions  of  0.01  ft/sec  at  low  velocities 
have  already  been  achieved  in  water. 

( d)  Direction  sensing  arrays  of  thermistors  can 
be  designed.  This  has  been  tried  in  air  and 

in  water  to  a very  limited  extent  (references  5* 
6,  and  7) • 

( e)  Changes  in  thermistor  characteristics  due  to 
changes  in  the  thermistor  itself  while  in 
operation  can  be  avoided  by  proper  seasoning  of 
the  thermistor  ( reference  l) . Changes  caused 
by  deposition  of  material  from  sea  water  on  the 
surface  of  the  thermistor  can  be  reduced  by  a 
flushing  or  cleaning  device  on  the  sensor. 

( f ) Physical  fragility  can  be  decreased  by  structural 
design. 

Therefore,  thermistors,  despite  development  diffi- 
culties, can  be  reliable  fluid  velocity  meters,  assuming  in 
water-  the  role  filled  by  the  hot  wire  in  air. 
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II.  BACKGROUND 

Davidson  Laboratory  has  developed  a thermistor 
velocity  meter  to  measure  bottom  velocities  under  ocean  swell 
(Figures  1 and  2),  which  has  been  used  successfully  to  measure 
fluctuating  velocities  up  to  0.2  ft/sec  peak  and  has  withstood 
several  days  of  intermittent  operation  in  water  depths  of  up 
to  forty  feet.  Deposits  on  the  thermistor  surface  during 
this  time  did  not  appear,  and  handling  the  probes  required 
only  those  precautions  normally  associated  with  a fragile 
instrument . 

The  voltage  output  from  the  thermistor  has  been 
amplified,  linearized,  and  written  out  using  the  circuit 
shown  in  Figure  Diode  linearizing  circuits  have  also  been 
designed  but  have  not  yet  been  tested.  Figure  3 shows  a 
typical  diode  circuit.  These  linearizing  circuits  have  been 
designed  (reference  5)  for  velocities  in  the  range  of  interest 
and  cam  be  made  to  linearize  an  input  as  accurately  as  desired 
by  dividing  the  thermistor  output  curve  into  smaller  and 
smaller  steps  and  by  providing  an  appropriate  diode  bias- 
voltage. 

Temperatures  of  the  fluid  near  the  velocity  sensor 
have  been  measured  with  a second  thermistor  (Figure  4),  When 
the  ambient  temperature  is  known,  the  correct  calibration 
curve  can  be  selected  fox-  the  reduction  of  the  velocity  data. 

Figure  5 shows  the  calibration  of  the  velocity  probe 
at  various  ambient  temperatures  for  our  oceanographic  work. 

The  thermistor  was  calibrated  by  towing  the  probe  at  different 
velocities  in  water  of  different  temperatures.  Figure  fc  shows 
a sample  of  oceanographic  data.  Figure  j0'  was  recorded  off 
Block  Island  In  August  19<ol  and  shows  the  simultaneous 
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recording  of  pressure,  temperature,  and  velocity  at  a point 
under  a wave  near  the  ocean  bottom. 

Figure  7 shows  a trace  taken  in  a wave  channel  to 
check  the  accuracy  of  the  thermistor  velocity-meter.  Surface 
elevation  was  recorded  as  a function  of  time  and  the  peak 
velocity  was  computed  using  shallow  water  wave  theory.  By 
examining  the  velocity  record  obtained,  the  peak  amplitude 
of  velocity  was  found  to  be  within  3 percent  of  the  computed 
value  of  approximately  0.5  ft/sec. 

Figure  9 shows  a trace  obtained  during  an  experi- 
ment in  a wave  channel  to  determine  the  phase  lag  (time 
constant)  of  the  thermistor.  A single  wave  of  short  period 
was  generated.  The  thermistor  and  a wave-height  indicator, 
which  used  a styrofoam  flat  linked  to  a rotary  differential 
transformer,  were  positioned  in  the  same  plane  perpendicular 
to  the  wave  front.  The  thermistor  was  in  the  potential  flow 
region  outside  the  bottom  boundary-layer.  By  comparing  the 
response  of  the  thermistor  to  that  of  the  wave-height  indi- 
cator, a close  approximation  of  the  time  constant  was  obtained, 
which  was  between  0.5  and  0.9  seconds. 

As  we  are  presently  using  them  to  measure  velocity, 
thermistors  are  insensitive  to  flow  direction.  However, 
arrays  could  be  designed  that  would  indicate  the  direction 
of  flow  ( reference  5) • Experiments  have  been  carried  out  at 
Davidson  Laboratory  using  a single  thermistor  as  a planar 
velocity  direction  finder  in  water.  Figure  8 shows  a data 
trace  obtained  in  the  wave  channel.  The  probe  was  aligned  with 
the  flow  and  provided  a varying  signal  as  it  was  rotated 
through  180  degrees.  In  a fixed  position,  a varying  signal 
was  also  obtained  corresponding  to  the  sense  of  the  velocity 
of  the  wave. 

The  possible  range  of  velocity  measurement  using 
thermistors  can  be  increased  by  using  larger  thermistors  with 
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larger  power  Input,  up  to  several  hundred  milliwatts,  or  by 
providing  circuits  that  accept  the  widely  varying  signal 
produced  by  a thermistor  over  a wide  velocity  range.  Recently, 
experiments  have  been  performed  that  showed  thermistor 
response  up  to  3 ft/sec,  and  their  response  probably  can  be 
extended  further. 
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FIGURE  2 


MECHANICAL  DRAWING  OF  THERMISTOR  VELOC I T Y -PROBE 


IGURE  4.  SCHEMATIC  DIAGRAM  OF  DIODE  LINEARIZING  CIRCUIT 
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FIGURE  5 


TEMPERATURE  BRIDGE  FOR  2000  OHM 
THERMISTOR  (FiNS  ARE  LETTERED  FOR 
SANBORN  INPUT) 


FIGURE  10  WAVE  CHANNEL  DATA  TO  DETERMINE  TIME  CONSTANT 
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FIGURE  n.  THEORETICAL  RESPONSE  OF  FENWAL  GB43p2  THERMISTOR 
WITH  RAT  ~ I06 
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